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Abstract. Star clusters are often used as tracers of major star formation events in external 
galaxies as they can be studied up to much larger distances than individual stars. It is vital to 
Q\ ■ understand their evolution if they are used to derive, for example, the star formation history 
of their host galaxy. More specifically, we want to know how cluster lifetimes depend on their 
environment and structural properties such as mass and radius. This review presents a theoretical 
overview of the early evolution of star clusters and the consequent long term survival chances. 
It is suggested that clusters forming with initial densities of > 10 4 Mq pc -3 survive the gas 
expulsion, or "infant mortality", phase. At ~ lOMyr they are bound and have densities of ~ 
00 ' 10 3±1 Mq pc -3 . After this time they are stable against expansion by stellar evolution, encounters 
' with giant molecular clouds and will most likely survive for another Hubble time if they are in a 
^v^j . moderate tidal field. Clusters with lower initial densities (< 100 Mq pc -3 ) will disperse into the 
field within a few 10s of Myrs. Some discussion is provided on how extra galactic star cluster 
populations and especially their age distributions can be used to gain insight in disruption. 
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i 1. Introduction 

5_i ■ In the Milky Way the majority of star formation occurs in embedded clusters (Lada 
& Lada 2003). At the same time, only a small fraction of stars in the disc resides in 
(open) clusters indicating that most of the clusters/associations are not long lived and 
disruption on time-scales of the order of the ages of open clusters is important. 
fNj ■ Historically, studies on the lifetimes of star clusters focussed on the open clusters in 
the Milky Way. In the end of the 50's the scarcity of old (> Gyr) open clusters was noted 
nearly simultaneously in three studies: van den Bergh (1957), von Hoerner (1958) and 
Oort (1958). It was already speculated then by these authors that this was due to their 
finite lifetimes. The first quantitative explanation for this came from Spitzer (1958) who 
showed that interactions with Giant Molecular Clouds (GMCs) are a plausible cause 
for the disruption of open clusters. In his model a cluster with a density of 1 M pc~ 3 
survives the periodic heating due to passing GMCs for about 200 Myr. Wielen (1971) 
constructed an age distribution of open clusters within 1 kpc of the sun and derived a 
median lifetime of 200 Myr. Because of this remarkable agreement with the theoretical 
predictions by Spitzer, GMC encounters have since then been held as the dominant 
disruption mechanism for open clusters. Another supporting argument for this idea came 
from the fact that the old open clusters are strongly concentrated towards the Galactic 
anti-centre, where the density of GMCs is low, and their survival chances high (van den 
Bergh & McClure 1980 and also Fig. EE). 

Later, Elson & Fall (1985) and Hodge (1987) constructed age distributions for the LMC 
and SMC clusters, respectively. They both compared to Wielen's result and showed that 
the age distributions in the Magellanic Clouds extend to older ages, and consequently 
concluded that these clusters survive longer than their Milky Way counterparts. 

This review attempts to summarise work done on the disruption of star clusters since 
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Figure 1. The distribution of young and old open clusters in the Galactic plane, based on the 
catalogue of Dias et al. (2002). The old open clusters 1 Gyr) are found preferentially towards 
the Galactic anti-centre. 



the seminal studies mentioned above. Three time spans are considered in which different 
destruction mechanisms are important: 

• first few Myr: Expulsion of residual gas from the star formation process 

• first few 100 Myr: Mass loss due to stellar evolution 

• first Gyr: Interactions with GMCs and 2-body relaxation in a tidal field 

These points will be discussed from a theoretical and observational viewpoint in § [2] and 
§[3j respectively. 



2. Theoretical considerations 

2.1. The first few Myrs: gas expulsion and "infant mortality" of clusters 

A lot of attention has recently gone to the disruption of star clusters just after their 
formation due to the removal of residual gas from the star formation process by the 
feedback from hot stars. This was already considered by Tutukov (1978) and Lada & 
Lada (2003) argue it could be the explanation for the strong drop at a few Myrs in 
the age distribution of Galactic embedded and open clusters. The arguments why this 
is destructive are as follows: if the embedded stars are in virial equilibrium with their 
surrounding, then the stellar velocities are too high once the gas is removed. The cluster 
consequently expands or even complete dissolves. If we define the total mass of the 
system as M to t = M g + M*, with M* the stellar mass and M g the gas mass, then virial 
equilibrium before gas expulsion can be written as 

2 GM tot 

v* a = . (2.1) 

2r v y ' 

Here G is the gravitational constant, t> rms is the root-mean square velocity of the stars 
and r v is the virial radius of the (embedded) cluster. The latter relates to the half-mass 
radius as r v w 1.25rh, slightly dependent on the density profile. Since v^ ms is three times 
the 1-dimensional velocity dispersion squared, af d , and r^ is 4/3 times the half-mass 
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radius in projection, r e g, Eg. I2.1l can also be written in observable quantities as 

2 GAftot 
<J id = > 

with r\ k, 10. In most models it is assumed that the gas expulsion time-scale, i cxp , is 
much shorter than the crossing time of the stars, i cr . The latter is defined as 

l - = 2V ~ 2 (g& ) (2 - 3) 

"°- 65Myr ( lO3At P c-0 ' (2 ' 4) 

For t CI > > t cxp the positions of the stars are static during the change of the potential 
and the cluster's response can be calculated using the impulsive approximation. 

In here the velocities of the stars right after gas removal are still defined by the total 
mass of the embedded cluster (Eq. 12. ip and hence the cluster will expand to find a new 
equilibrium. Lets take M* = eM to t, such that e is the star formation efficiency (SFE). 
The final radius can then be expressed in e and the initial radius, r v (0), as (e.g. Hills 
1980) 

T 6 

^fo) = 27^1- (2 ' 5) 

From this it can be seen that for e ^ 0.5 the cluster dissolves. Since this value is high 
for a SFE, gas expulsion is a plausible cause for the disruption of a large fraction of the 
embedded clusters. A value of e — 0.5, however, is for several reasons an upper limit for 
complete destruction of the cluster: 

• Escape of unbound stars: For e ^ 0.5 the system is globally unbound since the total 
energy (kinetic plus potential) is positive. However, most of the positive energy is carried 
away by escaping stars which have velocities much higher than u rms . Almost all dynamical 
models find that for e > 0.33 a fraction of the stars will remain bound (see Fig. [2]). 

• Sub-virial initial conditions: This was first discussed by Lada et al. (1984) who 
considered a collapsing cloud in which the stars are dynamically "cold" while the gas 
is removed (see also Proszkow et al. 2009). If the stellar velocities are only a fraction / 
of the virial velocities then Eq. 12.51 becomes r v /r v (0) — e/(2e — f 2 ), i.e. the condition 
for complete disruption becomes e = f 2 /2, which for / = 0.5 is e = 0.125 (Goodwin 
2008), much lower than the analytic estimate of e = 0.5 for initially virialised embedded 
clusters. The combination of e and the initial virial state, or the effective SFE (Goodwin 
& Bastian 2008), thus determines the survival chance. 

• Clumpy initial conditions: Fellhauer et al. (2009) showed that if stars form in a 
clumpy fashion, and if these clumps merge quickly, then the survival chance of the merger 
remnant is enhanced. 

• Slow gas removal: If the gas expulsion time-scale is longer than t cr , the cluster will 
expand adiabatically and this can result in a bound cluster for e as low as 0.1 (e.g. Geyer 
& Burkert 2001; Baumgardt & Kroupa 2007 and Fig. [2}. 

Therefore, the importance of the gas removal phase in disrupting star clusters depends 
on many parameters: the initial density of the embedded cluster, the gas removal time- 
scale, the initial virial state of the embedded cluster, the star formation efficiency and 
the initial configuration of the stars (clumpy or smooth density profile). This has led 
to several different predictions. For example: most models assume that gas removal is 
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Figure 2. The fraction of stellar mass that remains bound after gas expulsion for different SFEs 
as found by several authors using iV-body simulations. Data taken from Baumgardt & Kroupa 
(2007). 

instantaneous. If e is the same for all clusters then the disruption rate is independent of 
cluster mass. Baumgardt et al. (2008) advocate, however, that gas expulsion less efficient 
in disrupting massive clusters because of their higher binding energy per unit mass. This 
makes the general prediction that the exposed cluster mass function should be shallower 
than the embedded cluster mass function (see Parmentier et al. 2008). 

Since the embedded phase is short lasting (~ 1 Myr) most of these parameters are 
poorly constrained from observations and the interpretation of the models is therefore 
sensitive to the assumptions that are made for the initial conditions. In § [3] possible 
observational signatures of this early disruption will be discussed. 

2.2. The first few 100 Myrs: mass loss due to stellar evolution 

Clusters that survive the gas expulsion phase continuously lose mass through the evo- 
lution of the member stars. For a Kroupa type IMF between 0.1 M & and 100 M Q the 
cluster mass reduces roughly by 10%/20%/30% in the first 10/100/500 Myr. Since these 
stellar evolution time-scales are much longer than typical values of t cr , the cluster will in 
most cases react by an adiabatic expansion. This expansion is analogous to the increase 
of the orbital separation of a binary when one of the members loses mass. It is far less 
sensitive to the fractional mass loss than in the impulsive approximation, since 

r v (t) M(Q) 



r-v(O) M(t) 



1 

1*' 



(2.6) 
(2.7) 



In the last step \i is used for the fraction of mass that remains and this can be compared 
to e in the instantaneous mass loss case (Eq. 12. 5[) . For fj, = 0.5 the cluster will expand 
by a factor of two and not lose mass, whereas for impulsive mass loss the cluster nearly 
dissolves, as we have seen in § 12.11 

The situation is somewhat more complicated than this simple prediction. The expan- 
sion is more severe and can even result in complete disruption if the cluster is mass 
segregated before the bulk of the stellar evolution takes place (Vesperini et al. 2009). 
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If the cluster is not primordially mass segregated, it will still segregate while the stars 
evolve (Applegate 1986) and this can eventually lead to an enhanced expansion at later 
times. Finally, the fate of the cluster also depends on the Jacobi (or tidal) radius (rj). 
The combined effect of mass loss by stellar evolution and dynamical evolution in a tidal 
field was considered by Fukushige & Heggie (1995). They show that when clusters expand 
to a radius of ~ 0.5 rj they lose equilibrium and most of their stars overflow rj in a few 
crossing times. As was the case with gas expulsion, a lot more observational constraints 
on the initial conditions of star clusters, such as the degree of primordial mass segregation 
and the initial radius with respect to the tidal radius, are needed before we can make 
general statements about the effect of stellar evolution on the evolution of star clusters. 

2.3. Interaction with Giant Molecular Clouds 

An external disruptive factor, already considered by Spitzer (1958), is cluster encounters 
with Giant Molecular Clouds (GMCs). Since GMCs are typically more massive than 
clusters, the cluster is more affected by an encounter than the cloud (Theuns 1991). The 
cluster lifetime due to periodic heating of passing clouds is inversely proportional to the 
volume density of molecular gas, p gas , and proportional to the density of the cluster, a 
typical result for disruption by external tidal "shocks" (e.g. Ostriker et al. 1972) 

C « 1 Gyr ( P M ( Ph _ \ . (2.8) 

Here p® as w O.O3A/0pc~ 3 is the molecular gas density in the solar neighbourhood and 
the constant is taken from Gieles et al. (2006) and equals within a factor of two the 
original result by Spitzer (1958). The dependence of idis on p gas means that lifetimes of 
star clusters also roughly scale inversely with the observable surface density of molecular 
gas, Sgas enabling us to make order of magnitude predictions for the lifetimes of clusters 
in other galaxies. In spiral galaxies disruption by GMC encounters is especially important 
in the early stages of evolution since clusters form in the thin gaseous disc where p gas is 
high. Older clusters are typically more associated with the thick disc where p gas is low and 
GMC encounters are less frequent. Since young (<1 Gyr) clusters in spiral galaxies have a 
near constant radius (e.g. Larsen 2004), Gieles et al. (2006) argue that the dissolution time 
due to GMC encounters is longer for more massive clusters, due to their higher density. It 
is not clear whether this lack of a mass-radius relation is a universal property imprinted 
by formation, or the result of an evolutionary effect. This point will be addressed in §|3l 

2.4. 2-body relaxation in the tidal field of the host galaxy 

On somewhat longer time-scales stars are lost from clusters through 2-body relaxation in 
which stars get accelerated to the escape velocity, v csc (e.g. Ambartsumian 1938; Spitzer 
1940). The relevant time-scale of escape is, to first order, the half-mass relaxation time 
(Spitzer 1987) 

^1/2 3/2 

t lh = 0.138 h — , (2.9) 



VmG In A 

where N is the number of stars, m is the mean stellar mass and In A is the Coulomb 
logarithm: In A InO.lliV (Giersz & Heggie 1994). The relaxation time is approximately 
the time stars need to establish a Maxwellian velocity distribution. A fraction £ of the 
stars in the tail of the distribution have velocities larger than v CS c and consequently 
escape. Assuming that this high velocity tail is refilled every t T h then the dissolution 
time-scale is idis = irh/£e- For isolated clusters v csc = 2v rms . For a Maxwellian velocity 
distribution a fraction £ c = 0.0074 has v > 2v rms and then idis = 137 i r h- For tidally 
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limited cluster £ e is higher since v esc is lower. For a typical cluster density profile £ c » 
0.033, implying tdis ~ 30£ r h (Spitzer 1987). The escape fraction £ G is often taken constant 
(e.g. Gnedin & Ostriker 1997). But £ c depends on r h (through w rms ) and on the strength of 
the tidal field, or rj (through v csc ). Effectively, £ e depends on the ratio fh/fj ( e -g- Spitzer 
& Chevalier 1973; Wielen 1988). Gieles & Baumgardt (2008) show that £ e oc (r h /r T ) 3 / 2 
for rh/rj > 0.05 (the tidal regime). Together with Eg. 12.91 we then find for clusters on 
circular orbits in the tidal regime that idis oc N/oj, apart from the slowly varying Coulomb 
logarithm. Here u> = Vg/ Rg is the angular frequency in the galaxy and Rg and Vg are 
the galactocentric distance and the velocity around the galaxy centre, respectively. So 
for a flat rotation curve we find idis oc Rg for a cluster of a given mass (e.g. Chernoff & 
Weinberg 1990; Vesperini & Heggie 1997). This linear dependence of idis on Rg makes it 
difficult to explain the universality of the globular cluster mass function from dynamical 
evolution of a power-law initial cluster mass function (e.g. Vesperini et al. 2003), but this 

1/2 

will not be discussed here. Note that u> can also be written as u> oc p } , where pj is the 
density within rj, which is how Lee & Ostriker (1987) present it. 

Fukushigc & Heggie (2000) show that the linear scaling of tdis with N is slightly affected 
by the finite time it takes stars to find the Lagrangian "exit" points. Baumgardt (2001) 
showed that i cr is of importance and finds for equal mass clusters that t&i s oc t^tl^ 4 . 
Baumgardt & Makino (2003) find that this scaling also holds for models of clusters 
with a mass spectrum, stellar evolution and for different types of orbits in a logarithmic 
potential. Their result for idis can be summarised as 

(al , „ j „ yr (_J_) 3 ' 4 |£ (^g^) (1 - .), (2.10) 

where e is the eccentricity of the orbit. If the Coulomb logarithm is taken into account 
the scaling is something like idis oc TV 65 (Lamers et al. 2005). 



3. Observational signatures of cluster disruption 

3.1. The velocity dispersion of young clusters 

Ho & Filippenko (1996) were the first to demonstrate that one can empirically determine 
the dynamical mass, Md yn , of star clusters beyond the Local Group using their integrated 
light. They used high resolution spectroscopy to measure aid from the width of the 
spectral lines and high angular resolution (HST) imaging to measure r e s . The dynamical 
mass then follows from the virial relation (Eq. I2.2|) . They did this for the ~ lOMyr old 
super star cluster "A" in the dwarf starburst galaxy NGC 1569 and found that its Md yn 
is comparable to the masses of Milky Way globular clusters, concluding that cluster "A" 
is a young globular cluster. Since then this analysis has been applied to several dozen 
more clusters with different ages and masses and in different galaxies (e.g. Larsen et al. 
2004; Bastian et al. 2006; McCrady & Graham 2007; Mengel et al. 2008). 

3.2. Super-virial and dispersing? 

The empirically derived Md yn relies on the assumption that the cluster is in virial equi- 
librium. One way of verifying the validity if this assumption is by comparing Afdyn to the 
photometric mass, M p hot- The latter can be derived from the cluster's total luminosity 
and age and a comparison to an SSP model. Bastian et al. (2006) provide a compilation 
of Mphot and Md yn values of 19 clusters. They find that for the somewhat older clusters 
(> 50 — 100 Myr) there is good agreement between M p hot and Md yn - For many of the 
young (~ 10 Myr) star clusters Md yn >> Mphot- Goodwin & Bastian (2006) suggest that 



IAUS 266. Star cluster disruption 



7 




10 6 10 7 10 8 10 4 10 5 10 l 

T [y] M phot /r e(r [M e pcT 1 ] 



Figure 3. Left: Light to dynamical mass ratio for 24 clusters found in literature. The red and 
blue points refer to studies done in the optical and NIR studies, respectively. Different symbols 
correspond to different galaxy types. The lines show the photometric evolution of an SSP model 
in two filters. Right: Dynamical mass over photometric mass for the same clusters, shown as a 
function of M p hot/r e g, which is a proxy of aid fEq. !2.2[) . The full line is a prediction of the effect 
of binaries on Md yn (Eq. 13. f I) , with la and 2a variations due to stochastical fluctuations shown 
as dashed lines (Figure from Gieles et al. 2009). 

this is a signature of gas expulsion and the "infant mortality" scenario discussed in § 12.11 
i.e. they conclude that these clusters are super-virial. 

In the left panel of Fig.[3]an updated version of Fig. 5 in Bastian et al. (2006) is given. 
It shows the light to dynamical mass ratio of 24 young star clusters in different galaxies. 
The data are compiled from literature (McCrady et al. 2003 ; Bastian et al. 2006; Smith 
et al. 2006; McCrady & Graham 2007; Moll et al. 2007; Ostlin et al. 2007; Larsen et al. 
2008; Mengel et al. 2008) and will be discussed in more detail in Gieles et al. (2009) and 
Portegies Zwart et al. (2010). The observations taken in the visible are shown in blue 
and the ones taken in the NIR are shown in red. The photometric evolution in the V and 
K band from the Bruzual & Chariot (2003) SSP models using a Chabrier IMF is shown 
as a full (blue) line and a dashed (red) line, respectively. All clusters have an age close to 
lOMyr. This is because of observational constraints: younger clusters are more obscured 
and at lOMyr the red supergiants phase sets in making the cluster brighter (especially 
in the NIR, see Fig. [3]) and therefore easier to study. 

It is important to realise what initial conditions are required for the gas expulsion 
scenario to work. The models discussed in § l2.1l show that after instantaneous gas removal 
it takes roughly 10 initial crossing times to either completely disrupt or expand to a 
new virial equilibrium. In the models of Goodwin & Bastian (2006) this corresponds to 
~ 25 Myr since their clusters have i cr = 2.5 Myr in the embedded phase. Only with an 
initial t CI of a few Myrs it is possible to "catch" an unbound expanding cluster at 10 Myr. 
This corresponds to an initial density (stars and gas) of ~ 100M Q pc~ 3 (Eq. 12. 4[) . Since 
the clusters considered here have densities of ~ 10 3 M Q pc -3 (see Fig. |4j, gas expulsion 
can not be responsible for the high velocity dispersions. The densities in the embedded 
phase were at least a factor 1/e higher because of the gas, but probably even more 
since they have also expanded since then. A similar conclusion can be drawn from the 
measured velocities, which are of the order of 20kms~ 1 . With such velocities an unbound 



8 



Mark Giclcs 



cluster dissolves into the field in a few Myrs, incompatible with their radii being only a 
few pc. We conclude that the "infant mortality" phase of dense clusters only lasts a few 
Myrs. The clusters considered here have evolved for at least 10— 100 t cr and are therefore 
probably bound objects (an argument also made by McCrady and collaborators). 

3.3. Binaries 

The point remains that the measured velocity dispersions are too high for some of the 
clusters. It could be that the IMF assumption is invalid. The virial velocities are higher 
when the IMF is more bottom heavy, but the differences between the measured velocities 
and the virial velocities are probably too large to be explained by IMF variations. An 
alternative is that the velocities are dominated by orbital motions of binaries in the 
cluster, not taken into account in the virial relation (Eq. 12. 2p . Kouwenhoven & de Grijs 
(2008) considered this and they show that the Mdy n /Af p hot ra ti° can be elevated by 
several factors (see also de Grijs in this contribution). But they concluded that the 
contribution of binaries to the velocity dispersion is less then 5% for clusters with a%d > 
10 km s _1 . So binaries can not be responsible for the observed velocities according to their 
model. But Kouwenhoven & de Grijs (2008) did not include the mass dependent stellar 
mass-to-light ratio in their models and also did not take into account stars more massive 
than 20 M Q . So the effect will be stronger for the clusters considered in Fig. [3l since the 
light is dominated by massive stars (~ 15 M©) for which (primordial) multiplicity is high 
and short periods and mass ratio close to unity are common. 

Gieles et al. (2009) derive a simple analytic expression for the additional velocity 
dispersion due to orbital motions of stars in binaries, <7bin, relative to the dynamical 
dispersion of a virialised cluster, aid 



f_ 

.,, v 0.25 



(J_\i( mi W lQ 3 days \ f / M/r cS \ 1 
U.6/ \15M @ ) \ P J yiOHlQpc- 1 ) ' [ ' 

Here / is the binary fraction, q is the mass ratio, mi is the mass of the primary star (the 
contribution of secondary star is ignored) and P is the orbital period. The normalisation 
constants are based on properties of massive binaries in the Milky Way and appropriate 
for an age of lOMyr. Since the cluster light is dominated by the most massive stars 
the measured velocity dispersion is affected by binaries according to this relation. The 
prediction for -Md yn including the velocities of binaries (Afdyn oc a\ A + a^ in ) is shown 
in the right panel of Fig. [3] (full line). The dashed lines show the la and 2a spread 
in the model due to stochastics. The discrepancy between Md yu and M p h t can largely 
be explained by virialised clusters with a modest binary population. Assuming that all 
clusters host a similar binary population, the clusters with a low dynamical dispersion, 
for which we use M p h t/»"cff as a proxy (Eq. 12. 2[) . are affected most. This could be the 
reason why the data show a decreasing Mdyn/Afphot ratio for an increasing Mphot /'"eff 
ratio. 

This suggest that these clusters are survivors of the gas expulsion phase and are at 
the beginning of the stellar evolution phase. In the next Hubble time these clusters will 
adiabatically expand by a factor of ~ 2 due to mass loss by stellar evolution if they are 
not mass segregated (§ I2.2|) and it remains to conclude that these clusters are good young 
globular cluster candidates if the tidal field they are in is not too strong. 

3.4. The mass radius relation of star clusters 

Several studies have pointed out that there is a very weak dependence of cluster radius 
on mass/luminosity (e.g. Zepf et al. 1999; Larsen 2004; Bastian et al. 2005). For the 
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Figure 4. Mass-radius relation for the same 24 clusters shown in Fig. [3j using the same colour 
and symbol coding and r^ = (4/3)r c ff. Lines of constant half-mass density, = 3M/(8ivr^), 
are overplotted. The cluster densities imply 0.2 < i cr /Myr < 2 (Eg. I2.4|l . 



clusters in Fig. [4] this seems not to be the case. Although these are only 24 clusters, 
and it is far from a tight relation, these objects seem more consistent with a constant 
density of p h « 10 3±1 M 0) or r h w (M/10 4 M Q ) 1 / 3 within a factor of two. The difference 
between this figure and the studies mentioned above is that the clusters considered here 
are confined to a narrow (young) age range. The mass-radius relation is of importance for 
disruption since t^^ c oc ph f§ 12. 3p . So for a constant cluster density GMCs will destroy 
all clusters equally fast. If the mass-radius relation evolves in time, for example to a near 
constant radius as found by several authors, then t^f c will become mass dependent. In 
any case, the densities of these clusters imply that t^^f c is much longer than a Hubble 
time, even when p gas is 10 times higher than in the solar neighbourhood (Eq. 12. 8p . so 
GMC encounters will not play an important role in the disruption of these clusters at 
this stage. 

In a solar neighbourhood type tidal field the clusters in Fig.@]have r^/rj 0.03, within 
a factor of two true for all of them. This implies that these clusters are probably also 
very stable against the expansion due to mass loss by stellar evolution, since in § 12.21 it 
was argued that rapid dissolution due to mass loss by stellar evolution only occurs when 
r^/rj increases to ~ 0.5. A similar idea was recently put forward by Pfalzner (2009) who 
discovered two evolutionary sequences in young Galactic star clusters. A dense group 
that starts with a density of ~ 10 5 Mq pc -3 at a few Myr containing the Arches cluster, 
NGC 3603 and Trumpler 14, to which they refer as "starburst clusters". They seem 
to expand along a sequence of constant mass, where at 10 — 20 Myr there are the red 
supergiant clusters (RSGC01 and RSGC02) with a density of ~ 10 3 M pc -3 . The second 
sequence also starts at a few Myr, but with much lower densities (~ IOMqpc -3 ) and 
expands along a M oc 1/Vh sequence to densities comparable to the field star density. 
These are referred to as "leaky clusters". The clusters in Fig. 2] could be compared to the 
red supergiant clusters in the Milky Way i.e. the end stage of the starburst sequence. 
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Figure 5. Age distribution of SMC clusters based on the catalogue of Chiosi et al. (2006). The 
sample is split in small and the large clusters/associations, with the boundary at 6 pc. The 
histograms are made using a 0.5 dex bin width with different starting values (boxcar averaging). 



4. Discussion 

Up to the 90s the open clusters in the Milky Way and the populous clusters in the 
Magellanic Clouds were the targets of studies on disruption. HST has enabled us to 
establish the properties of larger populations of clusters containing more massive clusters 
in quiescent spiral galaxies (e.g. Larsen & Richtler 2000), interacting galaxies (e.g. Bastian 
et al. 2005; Whitmore et al. 1999) and merger remnants (e.g. Miller et al. 1997). 

The primary tool used for studies on the disruption of clusters is the age distribution. 
Different groups give varying weights to the various aspects described in § [2] in the 
interpretation of the results. There are two flavours of empirically based cluster disruption 
models: one is based on externally driven cluster disruption. In this model the dissolution 
time depends on the mass of the cluster and its environment. It was based on the age and 
mass distributions of luminosity limited cluster samples in different galaxies (Boutloukos 
& Lamers 2003). The variations in dissolution time-scales are explained by differences in 
the tidal field strength (Lamers et al. 2005) and the CMC density (Gieles et al. 2006). The 
second model assumes that internally driven cluster disruption is most important. In this 
model roughly 80-90% of the number of clusters is destroyed each age dex resulting in a 
(mass limited) age distribution that declines roughly as r _1 . In this model it is assumed 
that the infant mortality of clusters proceeds for a few 100 Myr and that the disruption 
rate is equally fast for all masses. This model is based on the cluster population of the 
Antennae galaxies (Fall et al. 2005; Whitmore et al. 2007) and it was argued by these 
authors that their model is universal and should be able to describe the age distributions 
of cluster populations in other galaxies as well. 

The universal cluster dissolution scenario has led to some controversy. Chandar et al. 
(2006) show that the age distribution of SMC clusters declines as r~ 85 , almost the 
same as in the Antennae. Gieles et al. (2007) used the same data set and show that the 
age distribution of massive clusters (> 10 3 ' 5 M Q ) is flat in the first few 100 Myr. They 
reproduce the r -0 85 result for the full sample and conclude that the decline is caused 
by detection incompleteness, confirmed by de Grijs & Goodwin (2008) using similar 
arguments. Boutloukos & Lamers (2003) show that for a constant formation history and 
no disruption, the age distribution of a luminosity limited sample declines as T~^ 1 ~ a ', 
where —a is the index of the cluster initial mass function. The index £ describes the 
cluster fading with age: F\(t) oc , with F\ being the flux at wavelength A of a 
cluster with constant mass. For the Z7(F)-band C ~ 1.0(0.7), resulting in t~ v , with 
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0-7 iSv ~ 1-0 due to fading alone when a = 2. If the cluster mass function is a Schechter 
function (Gieles 2009; Larsen 2009), the age distribution of a luminosity limited sample 
not affected by disruption is as steep as 0.9 < 77 < 1.4 (see also Konstantopoulos in this 
contribution). Both models rely on the assumption that the cluster formation history 
has been constant in the age range considered. Bastian et al. (2009) recently pointed out 
that this is invalid for interacting galaxies like the Antennae galaxies and that a recent 
burst can be confused by mass independent cluster disruption. 

Another issue, especially important for distant cluster populations, is what we actually 
define as a star cluster. Fig. [5] shows the age distribution of clusters and associations in 
the central region of the SMC (Chiosi et al. 2006). The sample is split in two sub-samples 
of ~ 200 clusters, based on their size. The age distribution of large clusters falls off much 
quicker than the one of the compact group. The median radius of both samples are 4.5 
pc and 9 pc, respectively. At a distance of 20 Mpc it will be very hard to tell these two 
groups apart. A speculative suggestion is that more of the short-lived "associations" or 
"leaky clusters" are included in cluster samples in distant galaxies, that would not be 
considered as genuine star clusters if they could be resolved. 
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